Power Law Distribution of Wealth in a Money-Based Model 
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A money-based model for the power law distribution (PLD) of wealth in an economically inter- 
acting population is introduced. The basic feature of our model is concentrating on the capital 
movements and avoiding the complexity of micro behaviors of individuals. It is proposed as an 
extension of the Equfluz and Zimmermann's (EZ) model for crowding and information transmission 
in financial markets. Still, we must emphasize that in EZ model the PLD without exponential 
correction is obtained only for a particular parameter, while our pattern will give it within a wide 
range. The Zipf exponent depends on the parameters in a nontrivial way and is exactly calculated 
in this paper. 
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I. INTRODUCTION 

Many real life distributions, including wealth alloca- 
tion in individuals, sizes of human settlements, website 
popularity, words ranked by frequency in a random cor- 
pus of text, observe the Zipf law. Empirical evidence of 
the Zipf distribution of wealth [1-9] has recently attracted 
a lot of interest of economists and physicists. To under- 
stand the micro mechanism of this challenging problem, 
various models have been proposed. One type of them is 
based on the so-called multiplicative random process [10- 
21]. In this approach, individual wealth is multiplica- 
tively updated by a random and independent factor. A 
very nice power law is given, however, this approach es- 
sentially does not contain interactions among individu- 
als, which are responsible for the economic structure and 
aggregate behavior. Another pattern takes into account 
the interaction between two individuals that results in a 
redistribution of their assets [22-25]. Unfortunately, some 
attempts only give Boltzmann-Gibbs distribution of as- 
sets [24,25], while some others [23], though exhibiting Zipf 
distributions, fail to provide a stationary state. 

In this paper, we shall introduce a new perspective to 
understand this problem. Our model is based on the 
following observations: (i) In order to minimize costs 
and maximize profits, two corporations/economic enti- 
ties may combine into one. This phenomenon occurs fre- 
quently in real economic world. Simply fixing attention 
on capital movements, we can equally say that two cap- 
itals combine into one.(ii) The disassociation of an eco- 
nomic entity into many small sections or individuals is 
also commonplace. The bankruptcy of a corporation, for 
instance, can be effectively classified into this category. 
Allocating a fraction of assets for the employee's salary, a 
company also serves as a good example for the fragmen- 
tation of capitals. Under some appropriate assumptions, 
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we shall establish a money-based model which is essen- 
tially an extension of the Eguiluz and Zimmermann's 
(EZ) model for crowding and information transmission 
in financial markets |26ll27j . The size of a cluster there is 
now identified as the wealth of an agent here. However, 
analytical results will show that our model is quite dif- 
ferent from EZ's [23] , which gives PLD with an exponen- 
tial cut-off that vanishes only for a particular parameter. 
Here, a Zipf distribution of wealth is obtained within 
a wide range of parameters, and surprisingly, without 
exponential correction. The Zipf exponent can be an- 
alytically calculated and is found to have a nontrivial 
dependence on our model parameters. 

This paper is organized as follows: In section 2, the 
model is described and the corresponding master equa- 
tion is provided directly. In section 3, we shall present 
our analytical calculation of the Zipf exponent. Next, 
we give numerical studies for the master equation, which 
are in excellent agreement with analytic results. In sec- 
tion 5, the relevance of our model to the real world are 
discussed. 



II. THE MODEL 

The money-based model contains N units of money, 
where N is fixed. Though in real economic environment 
the total wealth is quite possible to fluctuate, our as- 
sumption is not oversimplified but reasonable, given that 
the production and consumption processes are simulta- 
neous and the resource is finite. The N units of money 
are then allocated to M agents (or say, economic enti- 
ties), where M is changeable with the passage of time. 
For simplicity, we may choose the initial state containing 
just N agents, each with one unit of capital. The state 
of system is mainly described by n s , which denotes the 
number of agents with s units of money. The evolution of 
the system is under following rules: At each time step, a 
unit of money, instead of an agent, is selected at random. 
Notice that our model is much more concentrating on the 
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capital movement among agents rather than the agents 
themselves. With probability a/y/s, the agent who owns 
this unit of money is disassociated, here s is the amount 
of capitals owned by this agent and 7 is a constant which 
implies the relative magnitude of dissociative possibility 
at a macro level. After disassociation, this s units of 
money are redistributed to s new agents, each with just 
one unit. It must be illuminated that an real economic 
entity in most cases does not separate in such an equally 
minimal way. However, with a point of statistical view 
and considering analytical facility, this simplified hypoth- 
esis is acceptable for original study. Now, continue with 
our evolution rules. With probability a(l — 7/s), noth- 
ing is done. And with probability 1 — a, another unit 
of money is selected randomly from the wealth pool. If 
these two units are occupied by different agents, then 
the two agents with all their money combine into one; 
otherwise, nothing occurs. Thus, 1 — a in our model is 
a factor reflecting the possibility for incorporation at a 
macro level. 

One may find that as a is close to 1 and 7 is not too 
small, a financial oligarch is almost forbidden to emerge 
in the evolution of the system; but, if the initial state 
contains any figure such as Henry Ford or Bill Gates, 
he is preferentially protected. Note that the bankruptcy 
probability of moneybags is inverse proportional to their 
wealth ranks, and the possibility of being chosen is pro- 
portional to sn s /N; thus, the Doomsday of a tycoon 
comes with possibility an s j/N, which is extremely small 
for large s. Meanwhile, the vast majority, if initially poor, 
is perpetually in poverty, with no chance to raise the eco- 
nomic status any way. In addition, if middle class exists 
at first, it will not disappear or expand in the foreseeable 
future. Again, it may be interesting to argue that when 
a is slightly above zero, the merger process is prevailing 
and overwhelming, and all the capitals are inclined to 
converge. In this case, though the rich are preferentially 
protected, the trend in the long run is to annihilate them 
until the last. Of course, one-agent game is trivial. Like- 
wise, it is not appealing to observe the system when 7 
goes to and a to 1, since both merger and disintegration 
are nearly impossible-in other words, all the capitals are 
locked, thus the wealth pool is dead at any time. 

Following Refs.Hl EsL l2S| in the case of N > 1, we 
give the master equation for n s 
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has been used. We must point out that Eq.(l) is almost 
the same as the master equation derived in Ref.[27j for 
the EZ model except for an additional factor 7/s in the 
third term on the right hand side of Eq.(l). Notice that 
this term is significant because otherwise the frequency 
of the disintegration for large s agents would be too high. 

Now we introduce h s — sn s /N, which indicates the 
ratio of wealth occupied by agents in rank s to the total 
wealth, and a = /2(1 — a), that represents the maxi- 
mum ratio of the disintegration possibility to the merger 
probability in the whole economic environment. Then, 
one can give the equations for the stationary state in a 
terse form: 



and 



2(s + a) ^ 



h x = 



^h r h s 



(4) 



1 + a 



(5) 



According to the definition of h s , it should satisfy the 
normalization condition Eq.(3) 
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(6) 



When a is less than a critical value a c — 4 which will 
be determined numerically in section 4, one can show 
that Eqs.(4-5) does not satisfy the normalization condi- 
tion Eq.(3). This inconsistency implies that when a < a c 
the state with one agent who has all the N units of money 
becomes important |28l l29| . In this case, the finite-size ef- 
fect and the fluctuation effect become nontrivial and the 
master equations (1-3) is no longer applicable to describe 
the system [28|, |2£|. In this paper, we shall restrict our 
discussion to the case a > a r . 



III. ANALYTIC RESULTS 



When a > a c , one can show that h s 
ciently large s with 



V = 



A/s^ for suffi- 



(7) 



Notice that this equation is only consistent when rj > 2 
because otherwise the sum ^r=i r ^ r wou ld be divergent, 
and thus h s — > A/s v becomes an inconsistent formula. 

The derivation of Eq.(7) is described as follows: When 



3 



s is sufficiently large 

s s ~ 1 
h s — ■ - y h r h s — r 
2(s + a) ^ 

1 

_(£fc i _ r ^ + /i.O(- 5S _ T _)) 



TABLE I: The results of H for various value of a. 
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where 5 < 1 but is close to 1, 5{rj~ 1) > 1 and 2Si]—l—r] > 
1. Therefore 



which gives that as s — > oo 

/is = 

The value of J2^=i r ^r can be further evaluated: 
Introducing the generating function 



(9) 



G(X) = Xl ~ h r 



(10) 



r=l 



one can rewrite Eq.(4) as 
x(G' - /ii) + a(G - h lX ) = xG' + a(G - x) = xG'G 

or 

- 1) = a(G - x) (11) 

with the initial condition 

G(0)=0 (12) 

Since h s — > A/s'' as s -> oo, G is only defined in the 
interval \x\ < 1. From Eq.(6), we also have G(l) = 1. 
What we need to calculate is just 

oo 

G'(l) = ^V/i r 

r=l 

Since the left and the right hand sides of Eq.(ll) are both 
zero at x = 1, we differentiate both sides by x and obtain 

G"x(l -G) + G'(l - G) - xG' 2 = a(l - G 1 ) 

Let x — > 1 and one finds that G"(l — G) vanishes in this 
limit provided 77 > 2, thus 
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One immediately obtains that 

^r = g 

r— 1 

and the exponent 



1 - a/I -4/a 



(14) 



(15) 



G' 2 (l)-aG(l)+a = 



(13) 



which is a positive real number for a > 4. Notice that 
when a = 4, the exponent 77 = 2. This implies that our 
calculation is self-consistent, provided Eq.(6). In sum, we 
find from the master equation that h s obeys PLD when 
s is sufficiently large and a > 4. It may be important 
to point out that when s is small, h s also approximately 
obeys the PLD, and the restriction a > 4, introduced for 
the sake of discussing master equation, can be actually 
relaxed. This argument has been tested by the simulator 
investigation, which supplies the gap of analytical tools 
and verifies the analytical outcome. 



IV. NUMERICAL RESULTS 

We have numerically calculated the number 

oc 
r=l 

based on the recursion formula Eq.(4) with the initial 
condition Eq.(5). Table. 1 lists the results of H for vari- 
ous value of a. From Table. 1, one immediately find that 
the normalization condition is satisfied for a > a c = 4, 
which, again, indicates consistency of related equations. 

Fig. 1-2 show h s as a function of s in a log-log scale for 
a = 10, a = 4.5, respectively. From Fig.l, one can see 
that h s conforms to PLD for s » 1 with the exponent 
77 given by Eq.(15). Fig. 2 indicates that h s observes the 
Zipf law for nearly all s with 77 = 3.0. 
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FIG. 1: The dependence of h s on s in a log-log scale for 
a = 10. 
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FIG. 2: The dependence of h s on s in a log-log scale for 
a = 4.5. 



The fitted exponents for various values of a are plotted 
in Fig. 3. They are given by 

ln(/i900 / hiooo) 
ln(1000/900) 

Fig. 3 also exhibits the analytic results from Eq.(15). The 
analytic outcome fits the exponents calculated from re- 
cursion quite well for a > 4.2. However, when a — > 4.0, 
discrepancy is obvious, since the convergence of h s to the 
correct power law is then very slow. 

We have also performed computer simulation, which 
gives excellent agreement with theoretical results derived 
from Eqs.(4-5) for a = 8 and s < 10, see Fig. 4. For more 
about our simulator investigation and further analysis for 
a < 4, see Ref. [30|. 



FIG. 3: The calculated exponent r\ for different values of a. 
Black squares represent the numerical results of rj obtained 
from h a using the extropolation method, see text. The solid 
line represents the analytic result Eq.(15). 
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FIG. 4: h 3 for a — 8 from both numerical calculation and 
computer simulation. Black stars represent outcome of com- 
puter simulation for N = 2.5 x 10 s , 7 = 2 and a = 0.88889. 
Total 2 x 10 6 time steps were run and the final 5 x 10 5 time 
steps were used to count n s statistically. The circles represent 
the theoretical results derived from Eqs.(4-5). 



V. DISCUSSIONS 

In this paper, we have introduced a so-called money- 
based model to mimic and study the wealth allocation 
process. We find for a wide range of parameters, the 
wealth distribution n s ~ A/ s n+l with 77 given by Eq.(15) 
for sufficiently large s. The crucial difference between our 
model and the EZ model is that the dissociative proba- 
bility of an economic entity, after he/she is picked 
up, is proportional to 1/s in our model. However, the 
corresponding probability in the EZ model is simply pro- 
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portional to 1. This difference gives rise to divergent 
behaviors of n s . In the EZ model, n s ~ A/s 2 - 5 exp(— as) 
for large s [27j . When n s is interpreted as the number 
of individuals who own s units of assets, the choice of 
Td ~ 0(l/s) is reasonable. Actually, since at the first 
step, we randomly picked up a unit of money, the indi- 
vidual who owns s units of assets is picked up with a 
probability proportional to s. According to the obser- 
vation in real economic life, large companies or rich men 
are often much more robust than small or poor ones when 
confronting economic impact and fierce competition. If 
Td ~ 0(1), the overall dissociation frequency would be 
proportional to s which is totally unreasonable. 

In real economic environment, capitals and agents be- 
have similarly at some point. For instance, they both 
ceaselessly display integration and disintegration, driven 
by the motivation to maximize profits and efficiency. 
This mechanism updates the system every time, and 
gives rise to clusters and herd behaviors. Furthermore, in 
an agent-based model, it is usually indispensable to con- 
sider the individual diversity that is all too often hard to 
deal with. When it comes to the money-based model, 
this micro complexity may be considerably simplified. 
Finally, the conceptual movement and interaction among 



capitals is not as restricted by space and time as between 
agents. Therefore, when econophysics is much more in- 
terested in the behaviors of capitals than that of agents, 
it is recommendable to adopt such a money-based model. 

The methodology to fix our attention on the capital 
movements, instead of interactions among individuals, 
will bring a lot of facility for analysis; moreover, using 
such random variables as 7 and a to represent the macro 
level of the micro mechanism also help us find a possi- 
ble bridge between the evolution of the system and the 
protean behaviors of individuals. Whether the bridge is 
steady or not can only be tested by further investigation. 
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